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Abstract. We discuss the origin of X-Ray Flashes (XRFs), a recently discovered class of Gamma-Ray Bursts 
(GRBs). Using a simplified model for internal shocks we check if XRFs can be intrinsically soft due to some 
specific values of the parameters describing the relativistic outflow emerging from the central engine. We generate 
a large number of synthetic events and find that XRFs are obtained when the contrast rmax/Finin of the Lorentz 
factor distribution is small while the average Lorentz factor F is large. A few XRFs may be GRBs at large redshifts 
but we exclude this possibility for the bulk of the population. If outflows with a small contrast are commonly 
produced, even a large population of XRFs could be explained. If conversely the Lorentz factor distribution within 
the wind is broad, one should then rely on extrinsic causes, such as viewing angle effects or high redshift. 
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1. Introduction 

An intriguing discovery in recent years is the existence of 
a population of soft gamma-ray bursts (GRBs) with little 
or no emission above 50 keV fHeis e et al. 20011 Kippen et 
al. 2001. rSarraud et"'al. 2003) . These events, which have 
been called X-Ray Flashes (XRFs), share a number of 
characteristics with the classical GRBs (long duration, 
non-thermal spectra...), and there is now general con- 
sensus on the fact that XRFs represent an extension at 
low energies of the GRB population. In this context it 
is natural to verify whether the models developed to ex- 
plain the prompt emission of GRBs can also explain a 
population of soft bursts like the XRFs. This is a complex 
issue because XRFs can be explained either by extrinsic 
factors (e.g. viewing angle, redshift) or intrinsic factors 
(e.g. Lorentz factor, energy deposition...). An overview of 
the factors that could give rise to soft GRBs appeared 
in Zhang & Meszaros (2002). Recently several authors 
have discussed in detail the effects of the viewing an- 
gle on the softness of GRBs (e.g. lYamazaki et al. 20021 
[Zhang et al. 2004|ILamb et al. 2004|l . 

In this paper we concentrate on the impact of intrinsic 
parameters and we specifically address the following ques- 
tion: can the internal shocks model, which successfully ex- 
plains many properties of the GRB prompt emission, also 
explain XRFs without calling upon a particular set of ex- 
trinsic factors ? Our work is based on an analytical model 
that captures the essential physics of internal shocks. We 
demonstrate that internal shocks can produce XRFs quite 
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naturally and we discuss the conditions required for this 
to happen. 

The paper is organized as follows. Section |21 summa- 
rizes our current knowledge of XRFs. SectionOintroduces 
our analytical model of internal shocks. Section^jpresents 
the results of the simulation of a large number of GRBs 
and discusses the conditions required to produce XRFs. 
Section [31 summarizes our results and presents our conclu- 
sions. 



2. Observational properties of XRFs 

2.1. Gamma-ray bursts and X-ray flashes 

In 2001, Heise et al. reported the discovery of XRFs, 
short transients detected by the Wide Field Cameras of 
BeppoSAX in the range [2-26 keV] but not seen above 
40 keV by the GRBM on-board the same spacecraft (see 
IBoella et al. 1997l for a description of the BeppoSAX mis- 
sion). In order to clarify the relationship between XRFs 
and GRBs, Heise et al. (2001) compared the properties 
(duration and spectral hardness) of 9 XRFs with the X- 
ray counterparts of 16 GRBs also detected in the Wide 
Field Cameras of BeppoSAX. They concluded that "the 
statistical properties of XRFs display in all aspects a na- 
tural extension of the properties of GRBs" . Kipp en et al.[ 
(2001) analyzed the spectra of 9 XRFs simultaneously de- 
tected by BeppoSAX/WFC and BATSE and found that 
XRF spectra, like those of GRBs, are well fitted by the so- 
called Band function consisting of two smoothly connected 
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power laws IIBand et al. 1993)l . Defining x — E/Ep, where 
Ep is the peak of the E'^N{E) spectrum we have 

r a;"exp(-(2 + a)a;) if x < f^, 

N{E) cx < ^ / _fl\"-/3 (1) 
x' [ Iq:^ j exp(/3 — a) otherwise 

a and (5 being the two slopes respectively at low 
and high energy. The extended energy range of 
HETE-2 (from 2 keV to 400 keV) has allowed 
to further analyze the relationship between XRFs 
and GRBs (.BarrauTet al. 20031 IBarraud et al. 2004al 
ISakamoto et al. 2004[[Lamb et al. 2004|l . It appears from 
these studies that XRFs are a continuation of the class 
of long GRBs at low energy. The distribution of their 
duration is indeed consistent with long GRBs. Heise et 
al. (2001) showed that the tgo of 17 XRFs detected by 
BeppoSAX ranges from 10 s to 200 s and is compara- 
ble to the distribution of igo for the 36 GRBs studi- 
ed. The distribution of their spectral parameter a is 
also in agreement with the observed distribution for long 
GRBs. Kippen et al. (2001) showed that the distribu- 
tion of a in the sample of XRFs and GRBs they stu- 
died is consistent with what is found for bright BATSE 
bursts. Moreover, Barraud et al. (2004a) found values of 
a within the range predicted by the synchrotron shock 
models (—3/2 < a < —2/3), whatever the value of Eo 
{Eg being related to the peak energy Ep by the rela- 
tion Eo = Ep/(2 + a)). The distribution of the spectral 
parameter (3 is also comparable to the distribution of f3 
for long bursts with a mean value of —2.5. XRFs extend 
the well-known hardness- intensity correlation to soft, faint 
bursts ( |Kippen et al.| 2001, Barraud ct al. 2003). Lamb et 
al. (2004) and Sakamoto et al. (2004) showed that XRFs 
also follow and extend the Eiso - Ep relation discovered by 
Amati et al. (2002). It is therefore now generally accepted 
that the XRFs, XRRs (X-Ray Rich GRBs) and classical 
long GRBs form a continuum, and that they share a com- 
mon origin. 

Using GRBs detected by BATSE, Preece et al. (2000) 
found a distribution of the peak energy that is narrow and 
centered around 200 keV. With the discovery of XRFs, 
Heise et al. (2001) and Kippen et al. (2001) have shown 
that this distribution is broader than previously thought 
and that it is extended towards low energies, down to a 
few kcV. 

One of the first explanations proposed for the XRFs 
was that they could be GRBs observed at very high red- 
shifts (Heise et al. 2001). This hypothesis was however dis- 
carded by Barraud et al. (2003), in view of the similar du- 
ration distributions of XRFs and long GRBs. Additionally, 
the first upper limits and measured spectroscopic redshifts 
for XRFs contradicts the high redshift hypothesis with 
XRF 020903 at z = 0.25 ( |Soderberg et al. 2004| ), XRF 
040701 at z = 0.215 (Kelson et al. 2004), XRF 011030 
at z < 3.5 (Bloom et al. 2003), XRF 020427 at z < 2.3 
(Amati et al. 2004) and XRF 030723 at z < 2.1 (Fynbo 
et al. 2004). 



The remaining possibilities to explain XRFs are (z) 
GRBs with different intrinsic properties or (ii) standard 
GRBs viewed off-axis. In this paper we consider option 
(i) and study if the internal shock model of GRBs can 
also account for XRFs and we determine the conditions 
required to produce them. 

2.2. Defining an XRF 

While XRFs are best defined by their Ep, the photon e- 
nergy of the maximum of their vF^, spectrum, Ep is not 
always available for weak soft events. Consequently, we 
prefer to use the ratio i?x/7 of the 2 — 30 to the 30 — 400 
keV fluences to classify bursts into XRFs or GRBs. This 
fluence ratio is easier to compute and more robust than 
Ep, and it has been shown that it closely reflects the value 
of Ep when it can be measured l|Barraud et al. 2004ajl . 
Following Sakamoto et al. (2004) we consider XRFs those 
events with R^/^> 1. This definition calls for the following 
remarks. First, the separation between XRFs and GRBs 
is somewhat arbitrary since the present data do not show 
a bimodal distribution of -Bp (some authors call X-Ray 
Rich GRBs intermediate events with i?x/7 in the range 
0.3 to 1.). Second, the true fraction of XRFs depends on 
the definition of XRFs, but even more on the biases that 
affect their detection. Measuring the distances of a few 
XRFs could be a first step towards estimating their true 
fraction in a given volume. For instance, GRB 030329 (at 
a redshift of 0.168) would have been classified as an XRF 
at a redshift larger than z = 2. 

With the above definition, seven transients studied by 
Heise et al. (2001) and seven transients studied by Barraud 
et al. (2003) are XRFs. Some XRFs have been described 
in detail in the literature; they include XRF 020903 and 
XRF 030723 detected by HETE-2 (jSakamoto et al. 2004L 
Butler et al. 2004) and XRF 020427 detected by Beppo- 
SAX (Amati et al. 2004). Three ev ents, GRB 981226 
(jFrontera et al. 2000|l . GRB 990704 (jFeroci et al. 200111 
and GRB 000615 (Maiorano et al. 2004) classified as X- 
ray rich GRBs by the SAX team are XRFs according to 
our definition. 

3. XRFs in the context of the internal shock 
model 

3.1. A toy model for internal shocks 

The basic assumption of the internal shock scenario is 
that the central engine of GRBs is able to generate a re- 
lativistic wind with a highly non-uniform distribution of 
the Lorentz factor (with a contrast Fniax/rmin reaching 
at least a factor of 2). The observed emission is then pro- 
duced when layers of different velocities collide within the 
wind, the dissipated energy being radiated in the gamma- 
ray range by means of synchrotron shock emission (Rees 
& Meszaros, 1994). The evolution of this relativistic wind 
can be followed with a hydrodynamical simulation (Daigne 
& Mochkovitch, 2000) but this requires large amounts of 
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computing time which prevents one from considering a 
large number of cases and fully explore the parameter 
space. These detailed calculations have nevertheless shown 
that a simplified approach where the wind is represented 
by many shells which interact by direct collisions can also 
produce good results (Kobayashi et al. 1997; Daigne & 
Mochkovitch, 1998). The reason for this success is that 
the kinetic energy of the wind largely dominates over its 
internal energy so that pressure waves can be neglected in 
a first approximation. Going a step further we have deve- 
loped for this study a toy model where internal shocks are 
limited to the collision of only two shells of equal mass m. 
Obviously, we cannot obtain from this toy model any de- 
tailed information on the temporal profiles but we expect 
that the main features of the burst energetics will be pre- 
served. Shell 2 (of Lorentz factor is generated a time r 
after shell 1 (of Lorentz factor Fi < r2). This time inter- 
val (multiplied by 1 -|- z) represents an order of magnitude 
of the observed burst duration, since in the internal shock 
model time scales seen by the observer reflect the source 
variability. The average power injected by the central en- 
gine into the wind in this two shell approximation is given 

by 



E = 



mc 



(Ti + Ta) = Mtc^ 



(2) 



where M = 2m/T and f = (Fi + r2)/2 are the average 
mass loss rate and Lorentz factor. Shell 2 will catch up 
with shell 1 at the shock radius 



Ts = 2cr 



1^ 2 



p2 p2 
^ 2 ~ 1 



(3) 



The two shells merge at and the energy dissipated in 
the collision is given by 



-Bdiss = mc2(ri + r2-2rs) 



(4) 



where Fg = V^i^2 is the Lorentz factor of the shocked 
material in the merged shell. In order to produce a GRB 
this energy has to be radiated in the gamma-ray range 
with a characteristic broken power-law spectrum. If the 
synchrotron process is responsible for the emission, the 
peak energy (maximum of vF^) is 



Er, ^ Eg. 



(5) 



where B and Fe are the post shock magnetic field and 
electron Lorentz factor and Csyn = ^ • Assuming that 
a fraction of the dissipated energy is transferred to a 
fraction ( of the electrons we get 



(6) 



where ec^ is the dissipated energy per unit mass in the 
comoving frame. 

Similarly, if a fraction as of the energy goes into a 
disordered magnetic field generated behind the shock 



B = (STraspec^)^/^ 



(7) 



the peak energy can be written as 



Ep^CpTsp-ey 



(8) 



where p is the post shock density, x = 1/2, y = 5/2 and 

2 

(9) 



Cp = a,„(87rasc^)V2 (^^!!^ 

We have considered below the possibility that x and y 
can be different from 1/2 and 5/2 if for example the 
equipartition parameters arc not constant but vary with 
p and/or e. The possibility of non-constant equipartition 
parameters has been considered by Chevalier (2003) and 
used in afterglow modelling by Yost et al. (2003) who 
assumed that as varies with the shock Lorentz factor. 
For the prompt phase, Daigne & Mochkovitch (2003) 
have shown that the condition 2a; + y < 1 (which therefore 
excludes the standard values a; = 1/2 and y = 5/2) is 
often required to obtain good fits of the temporal and 
spectral evolution of GRB pulses. 

In our two shell approximation, the physical parame- 
ters of the shocked layer Vg, Fg, p and e can be related to 
the wind quantities E, t, T and the contrast of Lorentz 
factor K = F2/F1 in the following way 



8crF^ 



(k2-1)(k+1)2 



2F 



kV2 + K-1/2 



M 



E 



1 

l + - 

K 



(10) 



(11) 



(12) 



(13) 



Replacing Fg, p and e by their expressions in Eq.8 yields 



E"" (p^yjK) 
where the function 



(14) 



2x 



(At2 - 1) (1 + 1/k)^ («l/2+^-l/2_2)' 



k1/2 



-1/2 



(15) 



has been represented in Fig.l for k = 1 to 10 and three 
choices of x and y. 

In spite of the simplicity of the two shell approxima- 
tion, Eq.l4 predicts an anti-correlation between duration 
and hardness as observed in real bursts (Kouveliotou et al. 
1993). Another important (and surprising) consequence of 
Eq.l4 is that Ep is a decreasing function of F as long as 
X > 1/G. This can be understood from Eq.lO which shows 
that internal shocks occur closer to the source in a flow 
with a low Lorentz factor, due to a large baryon load. If 
X > 1/6 the reduced Lorentz factor cannot compensate for 
the resulting increase of p (Eq. 12) and a harder spectrum 
is produced. To obtain softer bursts, "clean fireballs" (i.e. 
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Fig. 1. The function (pxy(K) for x = 1/2 and y — 5/2 (full 
line) X = y = 1/2 (dotted line) and x = y = 1/4 (dashed 
line). In the first case, the large value of y leads to a very 
steep dependence of (pxy on k. 



with a large f ) are required. This however only applies to 
cases where T remains sufficiently high so that pair opaci- 
ty is unimportant at the source (Meszaros & Recs, 2000). 
With pair creation the situation becomes more compli- 
cated and has not been considered in this paper. 

3.2. Synthetic GRBs and XRFs 

The simplicity of the two shell approximation allows us 
to construct large samples of synthetic bursts to check 
if XRFs can be produced for some specific choice of the 
parameters. A synthetic event is determined by the wind 
dynamics which is fixed by the values of r, E, T and k, 
the spectral parameters Cp, x, y, a and /3 and the redshift 
z. For the spectral slopes we adopt a = —1 and (3 = 
—2.5 which correspond to the average values obtained in 
spectral fits of bright long GRBs performed by Preece et 
al. (2000) . We consider three different possible choices for 
X and y: {i) x = 1/2, y = 5/2, i.e. standard equipartion 
assumptions; {it) x = y = 1/2, if for example the fraction 
of accelerated electrons is proportional to e so that Fg 
remains approximately constant (Daigne & Mochkovitch, 
1998); {in) x = y = 1/A which was used by Daigne & 
Mochkovitch (2003) in their description of the temporal 
and spectral evolution of GRB pulses. These smaller va- 
lues of X and y would correspond to a situation where the 
dependence of the magnetic field or/and electron Lorentz 
factor on the dissipated energy is much weaker than in the 
standard case. We believe that such a possibility cannot be 



excluded in view of the uncertainties in the microphysics 
of the shocked material. 

We generated a large number of events (from a few 
thousands to one million) by making assumptions on the 
distributions of the burst parameters, basing on observa- 
tions or common hypothesis on the GRB physics and ori- 
gin. The best constrained parameters are the redshift and 
the duration. If long GRBs (and XRFs) are related to the 
explosive death of massive stars, their rate is directly pro- 
portional to the cosmic star formation rate ip* and their 
distribution in redshift can be deduced from V*(-2) for 
which we have adopted the analytical expression given by 
Porciani & Madau (2001) with a maximum at z ~ 1.5 
(their SFR 1). As shown by Bloom (2003) it remains 
presently impossible to decide between the three possible 
SFRs proposed by Porciani & Madau (2001) which can all 
be made compatible with the present GRB redshift data 
when corrected for high- redshift bias. 

The distribution of the observed duration Iqq for long 
BATSE bursts is approximately log-normal with a maxi- 
mum at fgo ^ 20 s. Wc have also adopted a log- normal 
distribution for r with a maximum at Tmax = 10 s and 
we checked a posteriori (see Sect. 4.3) that the resulting 
distribution of Tobs = (1 + z.)t for synthetic bursts agrees 
with that of fgo. 

The last four parameters T, k, E and Cp are much less 
constrained by observations and we simply take for them 
uniform distributions between 100 and 500 for F, and 1 
for Log K, 51 and 53.4 for Log E and Log Cp = Log C^°° ± 
0.5 where Cp^" is the value of Cp which produces a typical 
burst with Ep = 100 keV if £; = 10^^ erg.s'^, f = 300, 
K = 4, T = 5 s and z = 1. The upper limit of 53.4 for 
Log E has been estimated from the requirement that the 
synthetic LogiV-LogP relation agrees with the BATSE 
data (Stern et al. 2001). The comparison is shown in Fig. 2 
for 10^ synthetic events. 

The assumptions of uniform distributions for f , Log k, 
Logi^ and LogCp appear to be the simplest ones consi- 
dering our ignorance of the true distributions. The choice 
we have made also supposes that these quantities are inde- 
pendent, which may be wrong. Thus, we cannot expect to 
obtain from our results any reliable estimate of the pro- 
portion of XRFs relative to GRBs, but we can identify 
the range of wind paramaters that favors the production 
of XRFs. We will then know how the XRF/GRB ratio 
varies when the distribution of these parameters differs 
from our initial simplest choice. 

4. Results 

4.1. Ep distribution and softness - fluence relation 

We first obtained the Ep distribution of our synthetic 
bursts. The results are shown in Fig. 3 for the three con- 
sidered choices of x and y. The full line in Fig. 3 represents 
the distribution for the whole sample while the dotted and 
dashed lines respectively correspond to the sub-groups of 
bursts which would have been detected by BATSE and 
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Fig. 2. Differential peak flux distribution of BATSE 
bursts (Stern et al. 2001) compared to a synthetic dis- 
tribution (10^ events) with LogE = 53.4. 



HETE 2. A threshold of 0.2 ph s'^ cm-^ in the 50 - 
300 kcV energy range was assumed for BATSE while for 
HETE 2 we adopted 1 ph s"^ cm'^ both for FREGATE 
(between 30 and 400 keV) and the WXM (between 2 and 
10 kcV). These thresholds were estimated from the work 
of Band (2003) for the typical energy range of each ins- 
trument. They are only indicative and in practice also 
depend on the burst Ep and spectral indices. This effect 
has not been included in our analysis. It should be rather 
moderate for FREGATE and the WXM due to the rela- 
tive flatness of the sensitivity curves (Band, 2003). In the 
case of BATSE it will contribute to decrease the already 
small number of detected events at low energy. 

For X = 1/2 and y = 5/2 and x = y = 1/2 the dis- 
tribution of Ep for BATSE bursts is wider than the ob- 
served one, which is confined between 10 keV and 1 MeV 
(Preece et al. 2000). Conversely, the agreement is excel- 
lent for a; = y = 1/4 (since the value of Ep is then much 
less sensitive to the dispersion of the wind parameters). 
We have therefore adopted x = y = 1/4 in the remainder 
of this paper since it appears that this choice of x and y 
gives the best results both for individual bursts (Daigne 
and Mochkovitch, 2003) and statistically for a large popu- 
lation. 

Figure 3 shows that BATSE misses most of the low Ep 
events while HETE 2, which is less sensitive than BATSE 
in hard X-rays but has a lower energy threshold, can de- 
tect at least part of them down to ~ a few keV as was 
the case for Beppo-SAX. This is also illustrated in Fig. 4 
where we have represented the softness i?x/7 as a func- 
tion of the total (2 - 400 keV) fluence for a population of 
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Fig. 3. Distributions of Ep for bursts that can be de- 
tected respectively by BATSE (dotted line) and by HETE 
2 (dashed line). The whole sample is represented by a 
full line; top panel: x = 1/2, y = 5/2; middle panel: 
x = y = 1/2; bottom panel x = y = 1/ 4,. The curves 
are normalized so that the integral for the whole sample 

/p(iogi;p)diog£;p = i. 



1450 synthetic events which would have been detected by 
HETE 2. The total number of events produced was 3000, 
so that the detection fraction was about 1/2 (a smaller 
number of bursts was used in this case to avoid confusion 
in the figure). The two limits of the softness at 0.075 and 4 
respectively correspond to the hardest and softest bursts 
for which the two bands 2-30 and 30-400 keV are both in 
the low or high energy part of the spectrum. With the as- 
sumed values of a = —1 and (3 = —2.5 the softness limits 
are simply given by 



7 



]ldE_ 



= 0.0757 



for hard events and 
-^■"^dE 



-Rx/7 — 



r400 p_ 

ho E 



^■^dE 



3.957 



(16) 



(17) 



for soft events. The two horizontal dashed lines in Fig. 4 
separate the GRB, X-ray rich GRB and XRF domains. 
In agreement with Fig. 3 it can be seen that the model 
generates a population of XRFs that can be detected by 
HETE 2 (but would have mostly escaped detection by 
BATSE). Since, from the toy model, we know all the in- 
put parameters of these synthetic XRFs, we can identify 
the key ingredients necessary to produce them. For this 
purpose we compare below the distributions of z, r, E, V 
and K between XRFs and standard GRBs. 
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Fig. 4. Softness versus total fluence for a population of 
1450 synthetic events that can be detected by HETE 2. 
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Fig. 5. Rcdshift distribution of synthetic GRBs (dotted 
line) and XRFs (full line) obtained with our toy model. 

4.2. Redshift distribution 

Figure 5 shows that the redshift distributions of synthetic 
GRBs (dotted line) and XRFs (full line) are very similar 
(the adopted normalization is J P(log(l+^;)) dlog(l + 2:) = 
1). In the context of our simulation, XRFs are not stan- 
dard GRBs observed at large z. This was already strongly 
suggested by (i) the duration distribution of the observed 
XRFs which is comparable to that of long GRBs and (ii) 



the recent redshift determinations (or upper limits) ob- 
tained for several XRFs (see Sect. 2.1). Nevertheless, even 
at large redshift a bright GRB is still observable (Lamb 
& Reichart, 2000) but can appear as an XRF so that 
the XRF/GRB ratio is expected to increase with z. Fig. 5 
shows that this is indeed the case: the XRF/GRB ratio at 
z > 5 is more than 2 times larger than at z = 1. However, 
the number of events at large z is not sufficient to account 
for the bulk of the XRF population. We have checked that 
this was not a consequence of our specific choice for ip^, (z) 
which is maximum at 2; = 1.5. With SFR 2 of Porciani and 
Madau (2001) which is nearly constant at z > 2, the dis- 
tributions of both GRBs and XRFs remain similar, being 
only slightly shifted to larger z. In any case, most events 
stay confined between z = 1 and 5 and the majority of 
XRFs are not GRBs at large z. 

Our synthetic XRFs are therefore intrinsically soft due 
to some specific values of their relativistic wind parame- 
ters. 



4.3. Distribution of the wind parameters in GRBs and 
XRFs 

We compare in Fig. 6 the distribution of the four wind 
parameters (r, f , k, E) in synthetic GRBs and XRFs. 

Observed duration: 

The distribution of the observed duration t(1-|-2:) in GRBs 
and XRFs is represented in Fig. 6a. As expected, it is in 
good agreement with the BATSE duration distribution 
for long GRBs. The average duration of XRFs is approxi- 
mately 50% longer. This is not a consequence of a larger 
redshift but of a preferred longer intrinsic duration resul- 
ting from the duration-hardness relation Ep oc t~^/^ for 
X = 1/4 (Eq.l4). 

Average Lorentz factor: 

The average Lorentz factor in GRBs closely follows the 
uniform input distribution while large values of T are fa- 
vored in XRFs. This is again a consequence of Eq.l4 since 
Ep (X f-V2 for X = 1/4. 

Contrast of the Lorentz factor: 

The distribution of the contrast n shows a striking dif- 
ference between GRBs and XRFs. XRFs appears to be 
produced by relativistic winds where the contrast typi- 
cally docs not exceed a factor of 4. The maximum of the 
XRF distribution is located at k = 1.4. Conversely, the 
proportion of GRBs steadily increases with k. As a con- 
sequence of the small contrast of F in XRFs the efficiency 
for energy dissipation by internal shocks 



^1/2 ^ ^-1/2 _ 2 



mc2(ri -hFa) 



,1/2 



-1/2 



(18) 



is small, close to 1% at the maximum of the contrast dis- 
tribution. 

Injected power: 
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Fig. 6. Distributions of Tobs = (1 + z)t (a), T (b), Log k (c) and Log E (d) for GRBs (dotted line) and XRFs (full 
line). The thin dotted lines in (b, c, d) represent the uniform distributions used as input. The adopted normalizations 
are / P{Log Tobs) dLog Tobs = 1, /fgo" = 1; lo ^(^'^9 dLog k = 1 and J^f''^ P{Log E) dLog E = 1. 



The distribution of the injected power E shows that large 
E are favored in both GRBs and XRFs because events 
with low injected power often escape detection. XRFs are 
therefore not characterized by a deficit of injected power 
(even if the largest E are more frequently found in GRBs). 
They appear weak and soft due to the inefficiency of their 
internal shocks. The dissipated energy in XRFs and GRBs 



is compared in Fig. 7. The two distributions peak at 7 10^^ 
and 6 10^^ erg respectively. Again, this difference of nearly 
a factor of ten comes from the lower efficiency of internal 
shocks in XRFs relative to GRBs. Finally, the shock pa- 
rameters, through the value of Cp, also show some diffe- 
rences between GRBs and XRFs. A smaller Cp naturally 
favors the production of an XRF but a reduction of Cp 




Fig. 7. Distribution of the dissipated energy in synthetic 
GRBs (dotted hne) and XRFs (fuU hne). 

alone is not enough since, contrary to a low k, it increases 
the softness without simultaneously decreasing the radi- 
ated power (Barraud et al. 2004b). 

4.4. Amati relation 

Assuming that a constant fraction of the dissipated e- 
nergy is transferred to the electrons and radiated, it is pos- 
sible to check if our synthetic bursts follow the Amati rela- 
tion between the isotropic radiated energy and the value 
of Ep in the burst rest frame (Amati ct al. 2002). The 
results are shown in Fig. 8 for ag = 0.3, this rather large 
value of ag being required to maintain a reasonable overall 
efficiency 

/ = ae X /is (19) 

where /is is the efficiency of dissipation by internal shocks. 
Bursts which could be detected by HETE 2 have been 
represented by large dots in Fig. 8. For this sub-group, the 
best fit by a power law gives 

/ E \ ° 

= (to^ j 

The exponent is close to 0.5 as found in the observational 
Amati relation, which now extends over five orders of ma- 
gnitude in £^rad when XRF 020903 is included (Sakamoto 
et al. 2004). We however notice that the dispersion of syn- 
thetic bursts relative to the power law is larger than for 
observed bursts. If the small dispersion of the observa- 
tional Amati relation is confirmed in the future with a 
larger number of bursts it will provide a strong constraint, 
probably indicating that a new physical ingredient - such 



Fig. 8. The Amati relation for synthetic hurts. A sam- 
ple of 3000 synthetic events is represented, the large dots 
corresponding to those which can be detected by HETE 
2. The line is the best fit for the HETE sub-sample. 

as a correlation between some of the wind parameters - 
should be included in the models. 

5. Conclusion 

We have used a simple internal shock model to generate a 
large number of GRBs with different relativistic wind pa- 
rameters such as the average Lorentz factor f , the contrast 
K between the maximum and minimum Lorentz factor or 
the injected power E. We adopted a lognormal distribu- 
tion of the intrinsic duration and obtained the redshift dis- 
tribution assuming that the burst rate is proportional to 
the cosmic star formation rate. We also assumed standard 
values for the low and high energy slopes of the spectrum, 
a = —1 and (3 = —2.5 and discussed different possibilities 
regarding the shock parameters. Our aim was to identify 
the physical conditions leading to the formation of XRFs. 
We have found that our synthetic XRFs exhibit distri- 
butions of redshift, duration and injected power rather 
similar to those of GRBs but strongly differ in the distri- 
butions of r and K. XRFs are events where the contrast 
of Lorentz factors is small, predominantly between 1 and 
2. 

Since we do not know the true distributions of these 
parameters we cannot make any prediction on the relative 
fraction of XRFs and GRBs. With the uniform distribu- 
tions adopted for T and LogK we obtain 16% of XRFs, 
27% of X-ray rich GRBs and 57% of GRBs. This is not 
in agreement with the HETE results which show approxi- 
mately equal fractions of the three classes (Barraud et al. 
2004a). If however the weight of events with a low contrast 
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of Lorentz factors is increased, for example with a distri- 
bution of LogK which remains uniform but extends to 0.6 
only, so that the maximum contrast is 4 instead of 10, the 
agreement with the HETE results is improved with now 
31% of XRFs, 33% of X-ray rich GRBs and 36% of GRBs. 

If XRFs are produced by internal shocks with a low 
contrast in the Lorentz factor distribution, the efficiency of 
energy dissipation is expected to be smaller than in GRBs. 
This can in principle be tested from multiwavelength fits 
of the afterglow lightcurves and spectra. Such a study 
including for the first ti me an XRF (XRF 020903) was 
recently presented by Lloyd-Ronning & Zhang (2004) 
The efficiency for XRF 020903 indeed appears to be 
much smaller than what is found for the GRBs of 
the sample. However all the eflticiencies measured by 



Lloyd-Ronning & Zhang (2004) are quite large (0.1 for 
XRF 020903, between 0.4 and 1 for the GRBs) which can 
be very challenging for internal shock models. 

If conversely the Lorentz factor is always highly vari- 
able so that small values of k, rarely occur one should 
then look for another origin for XRFs such as viewing an- 
gle effects. This possibility has been studied in detail by 
Yamazaki et al. (2002, 2004). It supposes that the jets 
responsible for GRBs are uniform and have sharp edges. 
If these assumptions are verified, Yamazaki et al. (2002, 
2004) have shown that many of the observed and statisti- 
cal properties of XRFs can be accounted for by assuming 
a distribution of the jet opening angle such as c>c A0~^ 
(with 10"^ < M < 10^1 rad), XRFs being then obtained 
for large viewing angles 6^ > A9. 

Of course, the intrinsic and extrinsic origins for the 
softness of XRFs do not exclude each other. A low Ep can 
result from both a small contrast and a large viewing an- 
gle. Their respective contributions to the XRF population 
could be evaluated, for different distributions of k and A9, 
by adding the viewing angle effect to the present study. 

Acknowledgements. The authors would like to thank the 
anonymous referee for an extremely careful reading of the 
manuscript. His/her comments and remarks have strongly con- 
tributed improving our paper. 

References 

Amati, L., Frontera, F., Tavani, M., et al. 2002, A&A 390, 81 
Amati, L., Frontera, F., in't Zand, J.J.M., et al. 2004, A&A 

426, 415 
Band, D. 2003, ApJ, 588, 945 

Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 
Barraud, C, Olive, J-F., Lestrade, J. P., et al. 2003, A&A 400, 
1021 

Barraud, C, and the HETE-2 Team 2004a, Proceedings of 
the "2003 GRB Conference" held in Santa Fe, NM, USA, 
Sept.8-12, 2003 

Barraud, C, Daigne, F., Mochkovitch, R. & Atteia, J.L. 2004b, 
Proceedings of the "2003 GRB Conference" held in Santa 
Fe, NM, USA, Sept.8-12, 2003 

Bloom, J.S. 2003, AJ, 125, 2865 

Bloom, J.S., Fox, D., van Dokkum, P.O., et al. 2003, ApJ, 599, 
957 



Boella, C, Butler, R., Perola, C, et al. 1997, A&AS, 122, 299 
Butler, N., Dullighan, A., Ford, P., et al. 2004, 

astro-ph/0401020 
Chevalier, R.A. 2003, lAU CoUoq. 192, Supernovae: 10 Years 

of 1993J, eds. J.M. Marcaide, K.W.Weiler, Springer Verlag 
Daigne, F., & Mochkovitch, R. 1998, MNRAS, 296, 275 
Daigne, F., & Mochkovitch, R. 2000, A&A, 358, 1157 
Daigne, F., & Mochkovitch, R. 2003, MNRAS, 342, 587 
Feroci, M., Antonelh, L.A., Soffita, P., et al. 2001, A&A, 378, 

441 

Fynbo, J.P.U., SoUerman, J., Hjorth, J., et al. 2004, ApJ, 609, 
962 

Frontera, F., Antonelh, L.A., Amati, L., et al. 2000, ApJ, 540, 
697 

Heise, J., in't Zand, J., Kippen, R. M., & Woods, P. M. 

2001, Gamma-Ray Bursts in the Afterglow Era. Edited by 

Enrico Costa, Filippo Frontera, and Jens Hjorth. Berlin 

Heidelberg; Springer, 2001, p. 16 
Kelson, D.D., Koviak, K., Berger, E., & Fox, D.B. 2004, GCN 

2627 

Kippen, R. M., Woods, P. M., Heise, J., in't Zand, J., Preece, 
R. D., & Briggs, M. S. 2001, Gamma-Ray Bursts in the 
Afterglow Era. Edited by Enrico Costa, Filippo Frontera, 
and Jens Hjorth. Berlin Heidelberg: Springer, 2001, p. 22 
Kobayashi, S., Piran, T., & Sari, R. 1997, ApJ, 490. 92 
Kouveliotou, C, Meegam, C.A., Fishman, G.J., et al. 1993, 

ApJ, 413, LlOl 
Lamb, D.Q., & Reichart, D.E. 2000, ApJ, 536, 1 
Lamb, D.Q., Donaghy,T.Q., & Graziani, C. 2004, New Astr. 
Rev., 48, 459 

Lloyd-Ronning, N.M., & Zhang, B. 2004, ApJ, 601, 371 
Maiorano, E., Masetti, N., Palazzi, E., et al. 2004, Third Rome 

Workshop on Gamma-Ray Bursts in the Afterglow Era. 

Edited by Marco Feroci, Filippo Frontera, Nicola Masetti 

and Luigi Piro. ASP Conference Series, Vol. 312, p. 221 
Meszaros. P., & Rees, M.J. 2000, ApJ, 530, 292 
Porciani, C, & Madau, P. 2001, ApJ, 548, 522 
Preece, R.D., Briggs, M.S., Mallozi, R.S., et al. 2000, ApJS, 

126, 19 

Rees, M.J., & Meszaros, P. 1994, ApJ, 430, L93 
Ricker, G., Atteia, J-L., Crew, G., et al. 2001, AIP Conf. Proc, 
662, 3 

Sakamoto, T., Lamb, D.Q., Graziani, C, et al. 2004, ApJ, 602, 
875 

Soderberg, A.M., Kulkarni, S.R., Berger, E., et al. 2004, ApJ, 
606, 994 

Stern, B. E., Tikhomirova, Y., Kompaneets, D., Svensson, R., 

& Poutanen, J. 2001, ApJ, 563, 80 
Yamazaki, R., loka, K., & Nakamura, T. 2002, ApJ, 571, L31 
Yamazaki, R., loka, K., & Nakamura, T. 2004, ApJ, 607, L103 
Yost, S.A., Harrison, F.A., Sari, R., & Frail, D.A. 2003, ApJ, 

597, 459 

Zhang, B.& Meszaros, P. 2002, ApJ, 581, 1236 
Zhang, B., Dai, X., Lloyd-Ronning, N.M., & Meszaros, P. 2004, 
ApJ, 601, L119 



